Integration of the Friedmann equation for universes of arbitrary complexity. 
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An explicit and complete set of constants of the motion are constructed algorithmically for 
Friedmann-Lemaitre-Robertson- Walker (FLRW) models consisting of an arbitrary number of non- 
interacting species, each with a constant ratio of pressure to density. The inheritance of constants 
of the motion from simpler models as more species are added is stressed. It is then argued that 
all FLRW models admit a unique candidate for a gravitational epoch function - a function which 
gives a global time-orientation without reference to observers. The same relations that lead to the 
construction of constants of the motion allow an explicit evaluation of this function. In the simplest 
of all models, the ACDM model, it is shown that the epoch function exists for all models with A > 0, 
but for almost no models with A < 0. 



I. INTRODUCTION 

Recent cosmological observations |l| continue to pro- 
vide strong evidence that the ACDM model is a reason- 
able approximation to the recent history of our universe. 
However, the true nature of the dominant term - the 
dark energy - remains one of the foremost problems in 
all of science today With the emergence of at least 
three major players in current cosmology; the cosmolog- 
ical constant A (or, perhaps, a refinement), spatial cur- 
vature k and the matter M (dominated by the yet to be 
understood dark matter), the view that the history of 
our universe is best represented by a trajectory in phase 
space has become commonplace. Indeed, many (but no- 
tably not all [3]) analyses of the observations are now 
filtered through the ACDM model and reported by way 
of our current position in the f^A - plane. Here we use 
this phase space approach to completely solve the central 
problem in classical theoretical cosmology, the integra- 
tion of the Friedmann equation. This program is com- 
pleted for all models consisting of an arbitrary number of 
non-interacting species and is accomplished through the 
algorithmic construction of a complete set of constants of 
the motion. These constants are to be considered as the 
fundamental characteristics of the model universes. In 
this approach the universe is viewed as a whole without 
reference to any particular class of observer. The ques- 
tion then naturally arises as to whether or not events in a 
given universe can be ordered globally without reference 
to any particular class of observer. This question is also 
explored by way of the construction of a gravitational 
epoch function. It is shown that such a function need 
not exist, but all observations strongly suggest that we 
live in a universe for which it does. 



II. THE FRIEDMANN EQUATION 

We use the following notation: a{t) is the (dimen- 
sionless) scale factor and t the proper time of comoving 
streamlines. H = a /a, q = ~da/a^, ' = d/dt, k is the 
spatial curvature (scaled to ±1 for fc ^ and in units of 
(L signifies length)), p is the energy density and p 



the isotropic pressure. A is the (fixed) effective cosmolog- 
ical constant and o signifies current values. In currently 
popular notation, the Friedmann equation is given by 



where 



and 



3IP 



n = 



inGp 



c^k 



,2 ■ 



(1) 



(2) 



(3) 



Consider an arbitrary number of non-interacting sepa- 
rately conserved species so that 
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Pi 



with each species characterized by 

Pi = WiPi 



(4) 
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where Wi is a constant and distinct for each species. (In 
the real universe different species interact. The princi- 
pal assumption made here is that these interactions are 
not an important effect.) The conservation equation now 
gives us 
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where for each species Ci is a constant. We can formally 
incorporate A and k into the sum ^ by writing wa = — 1 
so that Ca = c^A/i and Wk = —1/3 so that Ck = —(?k 
and write ([!]) in the form 



(7) 



where the sum it is now over all "species" . Since Wi is 
assumed distinct for each species, there is no way to dis- 
tinguish, say. Lovelock Q (i.e. geometric) and vacuum 
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contributions to A. However, one can introduce an ar- 
bitrary number of separate species about the "phantom 
divide" at w = — 1. 

If Ci > for species i it follows that all the classical 
energy conditions hold [5j] for that species as long as 

-l<w,<l. (8) 
The situation is summarized for convenience in Fig[TJ 
A k M R 
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FIG. 1: The range in Wi is shown (in the box) for which 
the species satisfies the classical energy conditions. Here M 
stands for dust {wm = 0) and R for radiation {wr — 1/3). 



It follows from ^ that 

i 

and writing 1 + z = 6], it also follows that 

(l + z)3(l+-0f]„ 



E,- %o(l + Z)3(l+-.) 



(9) 
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As 1 + 2 ^ cx) (Big Bang or Big Crunch) it follows that 
Hi for all species except that with the largest Wi > 
-1. Then Q, ^ 1. As 1 + z ^ 0, f^, ^ for all 
species except that with the smallest Wi < —1. Then 
fii — > 1. From ^ and (fTU)) it follows that the evolution of 
the universes considered are governed by the autonomous 
(but non-linear) system Q 



where ' = -il + z)j^ 



= n,i{2q - 1) - 3w^) 



(11) 



III. SOLVING THE FRIEDMANN EQUATION 

We solve the system ([TT|) algorithmically by way of the 
construction of a complete set of constants of the motion. 
Consider all species or any subset thereof. To distinguish 
this latter possibility we replace ihy a. A constant of the 
motion, defined globally over the history of a universe, 
and not with respect to any particular observer, is defined 
here to be a function 



for which 



T =0. 



Consider the product 



(12) 
(13) 
(14) 



where the exponents Va are constants. Note that we do 
not make use of ([7]) at this stage. For to hold for 
(fT4| . it follows from ^ that 



(2q- i)y = 3W 



where 



V 



WaVa. 



(15) 



(16) 



If V" 7^ then according to (|T5|) q itself is a constant and 
so from (fTTjl for each species the evolution is given by 



l-2q+3Wi 



(17) 



We do not consider this case further here. If q is not 
constant then for T of the form p4|) to be a constant of 
the motion it follows from (flSl) that 



(18) 



V = 



and the constant itself then reduces to 



(19) 



Sums, differences, products and quotients of constants 
of the motion are of course constants as well and so the 
range in a must be examined. First consider two species, 
say a — a,b. It then follows form (jl8p that Wa = w\, 
which is not possible as we require that the species be 
distinguished by w. Next suppose that a ranges over 
more than three species. It now follows from that 
any such constant can be reduced to products and quo- 
tients of constants involving only three distinct species. 
As a result we need consider only those constants J- con- 
structed from three distinct species which we now do. 

The trivial case is that of a universe for which there are 
only three species, say a, 5, and c. Then T = fi^J" ri^*" rij^"^ 
and since we can always choose Va, (HH) reduces to two 
equations in two unknowns. The constant of motion is 
then determined up to the choice in Va- For example, 
with dust we have the species A, M and k (assumed ^ 
8], the ACDM model) and it follows immediately that 



(20) 



We can now substitute for flk from ([7]). 

More generally, consider n species, n > 3. If, say, 
species a and b are chosen to construct J- then there are 
n — 2 choices for the remaining species, the number of 
independent constants of the motion. Further, for ex- 
ample, there are (n — 2)(n — 3)/2 dependent constants 
of motion that follow which involve a but not b. In any 
event, as with the case of only three species, since we can 
always specify one exponent, (|18p always reduces to two 
equations in two unknowns. 
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As the foregoing algorithm makes clear, if species are 
added to a model, but none taken away, the constants of 
motion are inherited from the simpler model. For exam- 
ple, irrespective of what we add to the ACDM model, 
given by ((20)) . remains a constant of the motion. 

Now consider evolution in a three dimensional sub- 
space, say r2a, Jib, ric- We can construct the constant 
J'l = ^]a^^b''^^c^ and the constant T2 = fll'^fll''fld where 
d is any species other than a, b or c. Then with ([7]) we 
replace n^. All species other than a,b ot c that now 
enter are replace with the aide of constants constructed 
from that species and any two of a, b and c. In the sub- 
space then we have two independent relations that relate 
fia, J7c- The intersection of these surfaces in the sub- 
space defines the evolution trajectory of the associated 
universe in that subspace. 

To amplify the foregoing, consider, for example, dust 
and radiation. We have the species A, M, R and k (as- 
sumed ^ 0). Following the algorithm outlined above we 
immediately obtain the constants [13] 



^1 



M 



^4 - 



(21) 



(22) 



First note that the constant J-i is inherited from the 
ACDM model as discussed above. Moreover, only two of 
the above constants are independent. For example, J-3 — 
T\T/^ and T2 = T\T\. In, say, the k,M,R subspace 
the history of a universe is obtained by the trajectory 
defined by the intersection of J-^ with one of the other 
constants and with 17^ replaced by 1 — ilA in 
that constant [TTj . 

We end our discussion of constants of the motion here 
by noting that there are an infinite number of universes, 
of arbitrary complexity (but containing the species A, M 
and fc), for which 







(23) 



M 



throughout the entire history of the universe. The fact 
that we would appear to live in such a universe is a prob- 
lem for some cosmologists - the "flatness problem" . We 
take the view that this "problem" arises only when too 
few species are considered, in the limit a two species 
model with only M and k, in order to put the issue in 
perspective. 

Central to the discussion given above is the system 
(fTTj) . We now discuss how this system can be used to 
make a fundamental distinction between model universes. 



IV. THE GRAVITATIONAL EPOCH FUNCTION 



manifold. As stated in [l^], time-oriented is often weak- 
ened simply to time-orientable, and as pointed out in 
[l^ . the local time orientation is, with guesswork, ex- 
trapolated to the universe as a whole. Here we take a 
gravitational epoch function to be defined as a scalar field 
constructed from dimensionless ratios of invariants, each 
derivable from the Riemann tensor without differentia- 
tion, such that the function is monotone throughout the 
history of the universe. The purpose of an epoch func- 
tion is to allow the ordering of events without reference 
to any specific class of observers or, of course, coordi- 
nates. The existence of an epoch function, when unique, 
allows the rigorous time-orientation of an entire manifold 
without guesswork. This is why the existence of an epoch 
function is interesting. 

In a conformally flat four dimensional spacetime the 
maximum number of independent scalar invariants deriv- 
able from the Riemann tensor without differentiation is 
four. These are usually taken to be the Ricci scalar R 
at degree d = 1 and the Ricci invariants r{d — 1) for 
d = 2,3,4. (Degree (d) here means the number of ten- 
sors that go into the construction of the scalar. Let Rabcd 
signify the Riemann tensor, R^ the Ricci tensor, R the 
Ricci scalar (= R^) and the trace-free Ricci tensor 
(= R^ — ^S^). The Ricci invariants are deflned by (the 
coefficients are of no physical consequence as they derrive 
from the spinor forms of the invariants): r(l) = -^S^^S^, 
r(2) EE -^S'^S§S-, and r(3) ee ^S'^S^S^Sl The 
Kretschmann invariant, for example, is not used since 
^RabcdR"'""^ = R^ + 48r(l) in the conformally flat case.) 
The dimension of these invariants is L^^'^. In a Robert- 
son - Walker spacetime only two of these invariants 
are independent and in this spacetime it follows that 
r(l)2/r(3) = 12/7, r{lf/r{2f = 3, i?Vr(3) = 12/7? 
and (i?Vr(2))2 = 3/^ where 



^ r(l) 



(24) 



It is clear that / given by p4)) . or any monotone func- 
tion thereof, is the only epoch function that can be con- 
structed in a Robertson - Walker spacetime. Universes 
can then be fundamentally categorized as to whether or 
not they admit monotone /. 

For convenience consider T = //48. It follows directly 
that in any Robertson - Walker spacetime 



1, 



(25) 



Whereas T is deflned through a bounce, its representa- 
tion in the il notation is not. This is of no concern here 
since we are concerned with the behaviour of T strictly 
prior to any bounce should one exist. With the aide of 
PT|) . applied only to the species fc, it also follows that 



In rigorous texts on general relativity, for example [12 | 
and [l3[, spacetime is deflned (in part) as a time-oriented 



^ ^ {n.-q-ir ■ ^''^ 
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If we now apply the restriction ([5]) and use the full set of 
equations pT|) it follows that in addition to Q 



(27) 



Let us again consider the ACDM model fl^ . It follows 
from (USD that 



3n 



M 



and from ^ and (HZD that 



, r?M + 4»A . 



and finally from ((TT|) applied to M and A that 



(28) 



(29) 



(30) 



A/ 



As a result, T has a global minimum value of when 



(31) 



The question then is, what ACDM models intersect the 
locus ([31]) and therefore fail to have T monotone 15]? 

Since Qm > clearly all ACDM models with A > 
admit the epoch function Now consider A < 0. For 
k = condition §1^ holds at {flM^A) = (4/3,-1/3). 
For fc ^ use in ^ to define 



- 4)= 



(32) 



With k ~ ~1 choose a value of VIm in the range < 
51 Af < 4/3 and insert this value into ([5^ . The resultant 
•?^mt (which is negative) is the constant of motion associ- 
ated with the intersection of the associated integral curve 
with the locus ([3T|l with the intersection taking place at 
the chosen value of JIa/. For k = \ choose a value of 
in the range 4/3 < VIm < oo and follow the foregoing 
procedure {J-mt is now positive). However, with k — 1, 
as Hm — ^ oo the constant J-mt reaches the lower bound of 
16/27. There is then the range < T < 16/27 for which 
the associated integral curves do not intersect the locus 



((3T|) . We conclude then that for the ACDM models, all 
models with A > admit the epoch function (pS)) . but 
for A < only those with k = 1 and < < 16/27 do. 
There is no evidence that our universe lies in this latter 
category. Indeed, all the current evidence would suggest 
that we are very far away from it. 

We now make some general observations on the mono- 
tonicity of T, but restricted to the case fife = [l6]. In 
this case it follows from ((26|) that the monotonicity of 
T is closely associated with the monotonicity of q. For 
example, if we assume that q is monotone decreasing, 
then we can put a limit on the maximum w allowed for 
a monotone T. To do this note that for models with a 
Big Bang the initial value of q is {3w + l)/2 where w 
signifies the largest value of w (assumed > —1) for all 
species. It follows from (|26l) then that T can be mono- 
tone only for w < 1/3. In view of ^ this is a rather 
surprising restriction, but one which is in accord with all 
observations. 



V. SUMMARY 

A complete set of constants of the motion has been 
constructed for all FLRW models consisting of an arbi- 
trary number of separately conserved species, each with a 
constant ratio of pressure to density. These constants of 
the motion are to be considered as the fundamental char- 
acteristics of the model universes. The unique candidate 
for a gravitational epoch function has been constructed 
for all FLRW models. In the simplest of all models, the 
ACDM model, it has been shown that the epoch func- 
tion exists for all models with A > 0, for no models with 
A = 0, and for almost no models with A < 0. This 
function allows the global ordering of events without ref- 
erence to any particular class of observers or, of course, 
coordinates. 
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